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Abstract 

The problem of forward production of hadrons in heavy-ion colhsion at RHIC is 
revisited with modification of the theoretical treatment on the one hand and with the 
use of new data on the other. The basic formalism for hadronization remains the 
same as before, namely, recombination, but the details of momentum degradation and 
quark regeneration are improved. Recent data on the p/vr and p/p ratios are used to 
constrain the value of the degradation parameter. The pT spectrum of the average 
charged particles is well reproduced. A prediction on the px dependence of the p/p 
ratio at = 3.2 is made. 

PACS number: 25.75.Dw 
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1 Introduction 

Theoretical study of hadron production in the forward direction in heavy-ion collisions is 
a difficult problem for several reasons. The empirical fact that the proton-to-pion ratio is 
large at large rapidity implies that neither fragmentation nor hydrodynamics can be suc- 
cessful in describing the process of hadron production in that region. Recombination is the 
natural hadronization mechanism for large baryon/meson ratio, but the parton momentum 
distribution at low and large momentum fraction x (contributing to hadronic Feynman 
xp in the range 0.3 < xp < 1-0) in nuclear collision is hard to determine, especially when 
momentum degradation and soft-parton regeneration cannot be ignored. The use of data as 
input to constrain unknown parameters is unavoidable; however, that is also where further 
complexity arises. Data on forward production at ?7 = 3.2 ± 0.2 depend on both pi and 
Pti resulting in a smearing of the x distributions of the partons that makes phenomenol- 
ogy difficult due to the inter-connectedness of all aspects of the dynamical problem. The 
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problem was first studied in the framework of the recombination model in [H [2l [3] with the 
effects of parton regeneration taken into account [H [5]. The original data from PHOBOS 
show the distribution of charged particles at large rj, but without px measurement the value 
of xp cannot be determined |6i|. BRAHMS has measured both t] and pt dependences of 
charged hadrons [7j, but without particle identification the p/n ratio cannot be inferred. 
Very recently, there are preliminary data that indicate the p/vr ratio at r/ = 3.2 to be very 
large, ~ 4 at pt = 1-1 GeV/c, 0-10% centrality in Au-Au colhsions at = 62.4 GeV [8J, 
about 3 times higher than the prediction in [5]. The aim of this paper is to reexamine the 
problem of forward production and show that with appropriate changes in the treatment of 
degradation, regeneration and transverse momentum, there can be an understanding for the 
large p/ir ratio in the fragmentation region. 

In addition to the new data on p/ir ratio there is also a new presentation of the p/p ratio 
by BRAHMS for = 62.4 GeV, where the value of Rp/p ~ 0.02 is given [9|. That value 
differs from the value 0.05 inferred from the figure presented in [TD], which was the value used 
in [5]. The new values of Rp/Tr and Rp/p are consistent with the implication that there are 
more quarks or less antiquarks than what were obtained in [5]. That provides a hint for us 
to look for the area in the formalism where the treatment of degradation and regeneration 
may be improved. Regeneration is an effect that depends on momentum degradation in 
forward propagation, which in turn depends on the degradation parameter k that is not 
known except by fitting the data. With new data available, the whole procedure needs to 
be revised. In this paper we change the strategy of our phenomenology in order to take 
advantage of the additional constraints provided by the particle ratios. 

The formalism for forward production is basically the same as discussed in [U [5]. We 
describe its essence in Sec. 2, but with special emphasis on changes that are necessary 
to improve the treatment. In Sec. 3 momentum degradation and quark regeneration are 
investigated with significant changes from [U [S]- How the data on particle ratio can be used 
to constrain k is discussed in Sec. 4, followed by consideration of the transverse momentum 
in Sec. 5. Conclusion is given in Sec. 6. 

2 Basic Formalism for Forward Production 

In the recombination model (RM) [H [21 [3] hadron production can be described by the basic 
equations 

H^'^ix) = f^^^F^J^,{x„x„x,)Rpix„x„x„x), (1) 

^ J Xi X2 X3 

H^^{x) = [^^F^f{xuX,)R^{xuX,,x). (2) 

J Xi X2 

for proton and pion, respectively, where only one-dimensional consideration is needed for 
forward production, with x = Xp = Iphl ^ for hadron, and Xi being the momentum 
fractions of partons [U [5]. The recombination functions (RF), Rp and -R^r, depend on the 
wave functions of the hadrons and are summarized in [3]. The major task to render Eqs. 
([1]) and ([2]) useful is to determine the parton distributions and F^^ for the problem 
at hand. For forward production the largest contribution can be attained if the quarks 
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arise from different initial nucleons so that their momenta do not have to be shared among 
the quarks originating from the same nucleon. That means 0:2, Xs) depends on a 

factorizable product of independent quark distributions Fg'{xi) at momentum fraction Xi 
of an incident nucleon after z/j collisions with the target nucleus B; similarly, F^^{xi,X2) 
involves quark and antiquark distributions. If only 3 or 2 nucleons in the projectile A are 
considered in each collision, we can define the p and vr distributions from such sources as 
H!^^ and H^^^, which are then related to the overall distributions for AB collision by 

H^-ix, h) = j ^(3)B(,, 6, ,) . (3) 



Ht\-, b)=j 6, .) , (4) 

where h is the impact parameter. These formulas are derived in [1]. Clearly to describe 
Hf^^ and 

is a simpler problem than in Eqs. ([1]) and ([2]), since the corresponding 
parton distributions are for 3 and 2 nucleons, respectively, in the projectile. In Eqs. ([3]) 
and (HI) a denotes the inelastic cross section of nucleon-nucleon collision, and Ta{s) is the 
thickness function for a tube in A at impact parameter s. 

The recombination equation for the reduced projectile going through the target B is as 
in Eqs. and ^ 



^ J Xi X2 X3 



Hi^^^{x,b,s) = [ ——F^f^{xi,X2;\s-b\)R^{xi,X2,x) , (6) 

J Xi X2 

where F^^^^ (xi, X2, X3; \s — b\) is the 3-quark joint distribution after 3 nucleons transverse the 
target nucleus at impact parameter \s — b\ in B. We shall neglect the minor flavor dependence 
of nucleons and quarks in the following. Similarly, Fgf^{xi,X2; \s — b\) is the qq distribution 
after 2 nucleons go through B. In Ref. [4J F^^^^ is assumed to have a factorizable form. We 
now give a derivation of that form and in the process determine the appropriate average 
number of collisions. The same follows for Fgg . 

If each nucleon in the projectile nucleus A at s makes on average z/ collisions in B, where 

u ^ u^^ils- b\) = ^TB{\s-b\) 

1 — exp — (tTb(|s — 6|) 

then 3 nucleons make on average 3u collisions. Assuming a Poisson distribution in u, we 
have 

i/f^(x;6,.) = ^i/f^(x;i.)P3pM , (8) 
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where 



P.i^) = -^e-' . (9) 
Applying Eq. ([8]) to ([5]) the sum over u can be moved past the integrals and we can write 
F^^f{x,,X2,X3-As-b\)=Y.Fl^'J^''{x,,X2,x,-u)P3,{u) . (10) 

Now, u is the total number of wounded nucleons experienced by the target nucleus B, 
irrespective of how it is distributed among the incident nucleons. With three such nucleons 
that are independent, we have 

i^,?f(^l,^2,X3;Z.) = l Y: Vr7^r(^l)^;^(^2)F;^(^3) , (11) 

where the summation over z/j is constrained by J2i = ^i each starting from z/, = 0. Each 
term in the summand is a product of single-quark distributions in a proton that has under- 
gone Vi collisions with the target. They include the effects of degradation and regeneration 
to be discussed below. 

The use of in the Poisson distribution in Eq. ffTOl) is based on the assumption that all 
three nucleons in the projectile are lined up in the same tube at impact parameter \s — h\ in 
since otherwise the forward partons are not nearby in the transverse plane and are unlikely 
to recombine to form a proton. Thus the same v'''^ applies to each of the three nucleons. 
Substituting Eq. ffTTl) into ffTOl) and making use of the implicit 5y^y-^J^y.^+y.^ contained in the 
summation in ffTTl) . the sum over v can readily be carried out, yielding 

F(Jf(a:i,a;2,X3;|s-6|) = nF;(x.) , (12) 

1=1 

where 

oo 

F'^ix.) = E f:;%x.)pM) (13) 

with v being defined in Eq. ([7]) for a pB collision. Using Eq. (fT2|) in ([5]) and then in ([3]) we 
have reduced the proton production problem in AB collision to the only issue at hand, i.e., 
how the parton distribution F^{xi) is to be determined. 

For forward production we ignore the production of resonances and their decays. Proton 
is in the symmetric state in SU{2) x SU{2) for (spin, isospin). In 2x2x2 = 4-1- 2a-|-2s for 
qqq, the symmetric state (2^, 2^) -|- (2^, 2^) is 8 out of a total of 64 states, so the statistical 
factor in Rp{xi,X2,X3,x) is 1/8. For pion there is no change in Rt,{xi,X2,x) from that 
given in 

3 Momentum Degradation and Quark Regeneration 

The problem of forward production in pB collision has been treated in the framework of 
the valon model, which connects the bound-state problem of a static proton (in terms of 
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constituent quarks) with the structure problem of a proton in collision (in terms of partons) 
[21 O [11] . Without momentum degradation the quark distribution in a free proton is given 
by 

F,{x., Q') = £ dyG{y)K (^^, Q'^ , (14) 

where G{y) is the valon distribution, y being the momentum fraction of the valon, and 
K{z, Q^) is the quark distribution in a valon, both of which have been parameterized and 
updated in [12]. With momentum degradation in proton-nucleus collision both G{y) and 
K{z,Q'^) are modified, as described in [H |5|. However, we have come to the realization 
that Eq. (IT^ itself needs modification, a new development which we now describe from the 
beginning. 

A proton has three valons, which are the constituent quarks in the bound-state problem. 
When a proton wounds u nucleons in the target nucleus, it does not matter which of the 
the 3 valons causes the wounding; they can act independently. It is important to recognize 
the possibility that one of the valons may not undergo any momentum degradation, while 
the other two are responsible for causing u wounded nucleons in the target. Although the 
probability of that is low, the valence quark in the undegraded valon would have higher 
momentum. The point is that we should consider all possibilities, which can be expressed 
in the form 

where the dependence, shown explicitly in Eq. ( IT4l) . is suppressed because it is at some 
unspecified low value that is not of central importance here. G^{y') is the modified valon 
distribution due to degradation to be discussed below, together with the upper limit of 
integration. The Poissonian averaging of /i, the number of nucleons in B wounded by a 
valon, allows yU to be zero, while the total number of wounded nucleons is fixed at v. Thus 
the way that the valons are treated in a projectile nucleon is analogous to the way that the 
nucleons are treated in a projectile nucleus. 

If the momentum fraction that a valon retains after a collision with a nucleon in B is 
then after collisions the modified valon distribution is 

y'G',{y') = £ dyG{y)K'^S (^^ - k^^ , (16) 

which satisfies the normalization condition 

/ dy%{y') = J dyG{y) = 1 . (17) 



The solution of Eq. ( fT6i) is 



G'^iy') = K-^ G {^'''y' 
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It is clear that the maximum value of y' is because of the /z-fold degradation, thus setting 
the upper limit of integration in Eq. (ITSl) . Furthermore, the average momentum of the 
degraded valon is 

(y% = j dy'y'G'^{y') = K^{y) = \K^ , (19) 

where {y) is the average momentum fraction of a valon in a free proton and is 1/3. Thus Eq. 
f|T9|) expresses the effect of degradation in this simple model of multiplicative momentum 
loss of the sequential collision process. 

The valence quark distribution in a proton after v collision is as expressed in Eq. f|T5|) . 
but with K{z) replaced by the non-singlet component K^siz)-, which is specified in [12j. 
Due to the dependence of G'^{y') in Eqs. ffTSl) and ffT^ . the sum over fi in f[T^ acquires 
special significance at low /i, as remarked earlier before that equation. It is the ^ = term 
that renders the valence quark distribution at intermediate Xi insensitive to the value of k. 
In this respect our treatment here is an improvement over that in [ll[5]. 

For the regenerated sea quark distributions the earlier treatment can also be improved. 
In [S] the quark distribution K{z) in a valon is written in the two-component form 

K{z) = Kns{z) + L"{z) , (20) 

where L"{z) represents the regenerated sea quark distribution in a valon, including gluon 
conversion. The regenerated q distribution, F^{xi), for a nucleon making u collisions with 
the target is then as given in Eq. (1151) . but with K{z) replaced by L"{z). We now realize 
that such a convolution equation gives only a part of the total q distribution because the 
momentum lost by a nucleon after u collisions is not totally accounted for by that convolution 
equation. The average momentum loss of a nucleon as a fraction of the initial momentum 
after u collisions is 1 — (x)^, where 

We assume that the momentum loss is converted totally to uu + dd pairs without strange 
quarks. Thus the regenerated q {u or d) distribution for each nucleon in the projectile, F^{x), 
should satisfy the sum rule (with the subscript i on Xi suppressed) 

dxF^^ix) = 1 (1 - «:-) . (22) 
We adopt the approximate form for the x dependence 

F',{x) = Ml - xT , (23) 

so fy = {\ — K^)lA{n + 1). We shall use n = 7, since that is suggested by the q parton 
distribution of a free nucleon for x not too small and at low Q^. For the values of k and 
that we encounter below, Eq. ( l23i) gives values of F!^{x), for x > 0.2, far greater than 
those obtained by the convolution of Gl{y') with L"{x/y'), as determined in |5|; the latter 
is therefore neglected hereafter. 
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To summarize, for quark {u or d) distribution in pB collision after u wounded nucleons 
in S, we have 

f;(x) = f;(x) + f^(x) , (24) 

where Fg^{x) is the valence quark distribution given by Eq. ( |T5l) with Kns{z) in place of 
K{z), and F^{x) is the regenerated quark distribution given by Eq. (123!) . The antiquark 
distribution is, of course, just the second term in (l24l) . 

4 Particle Ratios 

Having obtained the modified quark distribution due to degradation and regeneration, we 
can now use Eq. fl2^ in f|T3l) for the ith nucleon, and then in f|T2|) for qqq distribution 
emerging from 3-nucleons colliding with target B at impact parameter Is* — h\. That result 
can then be used in Eqs. (I3l) and ([5]) to determine the x distribution of produced proton in 
AB collision. Exactly the same procedure can be followed to obtain the spectra of tt and p 
with appropriate use of the the q distribution for qq and qqq recombination. 

We show in Fig. 1 the results of our calculation of the x distributions of p, vr, and p for 
K = 0.7, and b = 3.3 fm for 0-10% centrality in Au-Au collisions. The value of k is chosen for 
reasons to be given below. Evidently, the p distribution is much higher than the other two 
for X > 0.5, since it is due to the recombination of three valence quarks from three different 
nucleons in the projectile A. Moreover, it decreases more slowly with increasing x due to 
the slower decrease of valence quark distribution compared to the sea quarks. Thus the p/n 
ratio is large and increases with increasing x. The tt distribution is much higher than the p 
distribution, because of the effect of valence quark in tt that is lacking in p. Similar plots 
can also be made for other values of k, but in the absence of any data on the hadronic x 
distributions the comparison among different k values can better be presented in a different 
format, as shown below. The general trend is that lower k leads to higher level of q and 
therefore higher vr and p at low x. 

Recently, data have become available on the particle ratios of both p/p [9\ and p/vr 
It is then very revealing for us to make parameteric plots of those ratios for various values of 
X and K. We use Eqs. ([T]) and to calculate H^^{x, k) for b = 3.3 fm and h = p,'7i,p and 
show their ratios Hp/ Hp versus Hp/H„ in Fig. 2, in which the grid lines are for constant x 
(in solid lines) and constant k (in dashed lines). It is clear that all lines have negative slopes 
in that figure because q is involved in the numerator of Hp/ Hp, but in the denominator of 
Hp/H-,^. Large values of Hp/H^^ can be achieved only when x is > 0.4; that is the region 
where the valence quarks dominate and the sea quarks are suppressed. At fixed x both 
ratios depend sensitively on k, more so for Hp/ Hp than for Hp/H„, because of the number 
of q involved. The smaller k is, the more degradation there is, and the regenerated q boosts 
Hp/ Hp and suppresses Hp/H^,. 

The data on Rpip and Rp/.,, depend on the values of pt at which the hadrons are included 
in the determination of the ratios. Rhi /h cannot be identified with H^i / H^ until after the pt 
distribution is considered, a topic to be discussed in the next section. So far we have only 
treated the dynamical processes that lead to the x distributions. At fixed rj the longitudinal 
Pl and the transverse pt are, of course, not kinematically independent. The range of x that 
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is phenomenologically relevant to our study should correspond to the range of pt in which 
the experimental values of the particle ratios are determined. Since the mismatch between 
Rh'/h and Hh' / Hh is not large, as we shall discuss later, let us here mark on Fig. 2 the data 
point that corresponds to |8j, I9j 

Rp/^ = 4.08 ± 0.2, r/ = 3.2 ± 0.2, 0.9 < pr < 1-3 GeV/c , (25) 



i?p/p = 0.0231 ± 0.0012, y = 3.0 ±0.1, 0.5 < < 1-4 GeV/c . (26) 

The grid lines in Fig. 2 then suggest that the relevant values of x and k, are 

X ~ 0.55 and k ^ 0.67 , (27) 

For that reason the x distributions in Fig. 1 are shown for k = 0.7. 

What we have done so far is essentially the first step of an iteration process, in which the 
focus is on the x distribution. The next step is to consider the transverse momentum based 
on the result of the first step and to improve on the overall phenomenology. 



5 Transverse Momentum 

The Pt dependence of the produced particles has been discussed in [S]. Let us first give 
a summary of that consideration. Since hard scattering is suppressed in the fragmentation 
region, we ignore shower partons for x > 0.2. This approximation is supported by the data on 
the Pt distribution of charged particles at rj = 3.2 [7j, which shows an exponential behavior 
for Pt up to 2 GeV/c without up-bending due to power-law behavior. Thus we write the x 
and Pt distributions of a produced hadron h in the factorizable form 

' H^ix,K)V,{pT), (28) 



Pt dxdpT 

where the transverse part is normalized by 

' dpT PT VhipT) = 1 , (29) 



rendering 



= Hh{x, k) , (30) 



which is our starting point in Eqs. ([T]) and ([2]). 

The properties of Vh{pT) described in [5J are adapted from the treatment of pt distri- 
bution in central collisions at mid-rapidity for which the only recombination process is in 
the transverse plane [13]. Here, we have treated in detail the degradation, regeneration and 
recombination of partons in the forward production, so it is inappropriate to append a sep- 
arate recombination of thermal partons with independent recombination functions for the 
transverse component. Since no shower partons are involved, we shall simply take a common 
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exponential form for all liadrons, but allowing the inverse slopes Th to differ for hadrons with 
different masses, as suggested by hydrodynamical flow. Thus we write 



VhM = ^f-'^^"^' (31) 
with normalization chosen to satisfy Eq. fl2^ . We parametrize by 

n = To + mh{vt)\ (32) 

where the second term expresses the flow contribution. Since at large x the dominant momen- 
tum direction is longitudinal, the mass- dependent component of the transverse momentum 
is expected to be small compared to the thermal component characterized by Tq. 

Although X and pt appear independent in Eq. (l28!l . they are kinematically constrained 
when 7] is fixed. They are related by 

X = sinh T] . (33) 

At f] = 3.2, the range 1 < Pt < 1-5 GeV/c corresponds to 0.39 < x < 0.59. On the other 
hand, if the rapidity y is fixed, the relationship depends on the particle mass. At y = 3.0, 
the range 1 < Pt < 1-5 GeV/c corresponds to 0.32 < < 0.48 and 0.44 < Xp < 0.57. The 
value X ~ 0.55 determined from our theoretical grid lines in Fig. 2 lies within the range of 
X values above for the data on Rp/^ a.t t] = 3.2 in Eq. fl25l) and also within the range of Xp 
values above for the data on Rp/p at y = 3.0 in ( l26l) . This is a non-trivial achievement, since 
the formalism described in Sec. 3 makes no reference to pr, so the grid lines for the ratios of 
Hh'/h{x, k) at constant k, and x need not imply any px values that correspond to the relevant 
X and Pt values of the experimental Rh'/h at fixed rj. 

The Pt distribution given in [7] is to be identified with our calculation as follows 

^Hf,{x)VH{pT) , (34) 



27ipTdpTdr] 271 

since upon integration over pTclpTclcj) and using Eq. fl29|) it yields Hh{x). Strictly speaking, 
holding rj fixed on the LHS is not the same as holding x fixed on the RHS. But the data 
are analyzed at t] = 3.2±0.2, so there are bands of rj and x values in which Eq. flMl) is 
approximately valid. The data on Rp/nipr) are then to be related to our calculation by 



Hp[x{pT)] VpjpT) 

HMpt)] K(pt) 



Rp/Apt) = 7l:r(! :rr. , (35) 



where the ratio Hp/Hj^ is to be determined by fixing rj = 3.2 and k = 0.67. In Fig. 3 we 
show that ratio by the dashed line, which has a significant pt dependence. Furthermore, the 
average magnitude of Hp/Ht, accounts for the major part of Rp/-,,, and cannot arise without 
a realistic treatment of degradation and regeneration. The reason for the dashed line to 
increase with pt is that at fixed t] higher pt means higher x, where q is suppressed compared 
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to g, resulting in vr being suppressed relative to p. The solid line in Fig. 3 includes the effect 
of V^/VJr, which we get from Eqs. (13T|) and (!32|) 



VpjPT) 




exp 




(36) 



Since the mh{vt)'^ term in Eq. (!32|) is small compared to Tq, as shall show presently, the above 
ratio is approximately exp[(mp — m^)pr(f^t)^/7'o], which shows the effect of mass difference 
in elevating the dashed line to the solid line. The result of fitting the data [8jj gives 

{vtY/T^ = 0.7 (GeV/c)-2 . (37) 

This is consistent with mp{vtY Tq, when Tq is 0.2 GeV/c to be determined below. 

The pt distribution itself is an additional test of our model, since the absolute normal- 
ization is not canceled as in a ratio. The data [7] are for all charged hadrons without particle 
identification, for which we treat {h^ + h~)/2 as = [p + p + 1.2 (vr^ + tt")] /2, where 
the K/tt ratio of ~ 0.2 is used [9]. As the third step in our iteration process, we calculate 
Hh±{x)Vh±{PT)/'2TT holding x and k fixed as in Eq. ( 1271) and adjust Tq to fit the data according 
to Eq. (jUj). The result is shown in Fig. 4 for Tq = 200 MeV; it agrees with the data very 
well. Since the normalization is fixed by the Hh{x, k) functions and is not adjustable, a good 
fit is remarkable. 

Putting the obtained value of Tq in Eq. (!37|) . we have 

To = 0.2 GeV/c, {vtf = 0.028. (38) 

The value of (vt) — 0.17 seems reasonable in view of the dominance of longitudinal expansion 
in the fragmentation region. The significance of this work is, of course, not in the trans- 
verse aspect of the problem, but on the longitudinal momentum distributions of the forward 
particles, which affect the pt distribution. The large p/n ratio found in the BRAHMS data 
at T) = 3.2 cannot be understood without a proper treatment of the x distributions in the 
fragmentation region. 

As a prediction of this work, we can calculate the pt dependence of the p/p ratio at fixed 
1] = 3.2. Since Vp/Vp = 1 for /i = p,p, only Hh[x{pT)] contributes to the ratio Rp/p{pT)- 
Using Eq. (!33|) . we have 

HpjxjpT)] . . 

- hmM ■ ^^^^ 

The result for ^/s = 62.4 GeV, b = 3.3 fm and k = 0.67 is shown in Fig. 5. The range of pt 
covered by the plot corresponds to x roughly between 0.3 and 0.6 at rj = 3.2. Note that the 
result is for fixed r], not fixed y. The reason for the decrease of Rp/p{pT) with increasing p^ 
is clearly the increase of x where q at momentum fraction Xi, approximately x/3, becomes 
more suppressed than q at the same Xi. A verification of this prediction would lend further 
support to our model. 
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6 Conclusion 



This work differs from the earher attempt in ^ in three important ways. Firstly, new data are 
available that put more stringent constraints on unknown parameters. Secondly, significant 
modifications have been made in the treatment of degradation, regeneration, and transverse 
momenta. Thirdly, the order that phenomenology is carried out is reversed due to the new 
empirical knowledge about the particle ratios. Using p/vr and p/p ratios as input, we are able 
to determine the degradation parameter k, which enables us to calculate the x distributions 
of the hadrons. At fixed 7] that implies a pt dependence of the p/vr ratio arising from the 
X distributions; that pt dependence accounts for a large part of the data on that ratio, the 
balance being due to the exponential pt distributions that are mass dependent. In fitting the 
particle ratio the calculated result is insensitive to the absolute normalization of the yield. 
The latter is shown to be correct when we succeed in reproducing the pt spectrum of the 
average charged particle. That is a significant achievement because the yields of protons, 
pions and antiprotons at large r] depend strongly on the dynamical process of momentum 
degradation and soft-parton regeneration. 

Although the degradation parameter k is determined by data fitting, to get the spectra 
correctly for all hadrons through one such parameter relies on the validity of the treatment 
of the various subprocesses. Our results suggest that our model has captured the essence of 
the dynamics involved. In particular, the large p/vr ratio would not have emerged from our 
calculation if recombination has not been used as the mechanism for hadronization. 

Since proton production at large x is due to the recombination of three valence quarks 
from three nucleons in the projectile, in which there are numerous other valence quarks 
from other nucleons, we do not expect the events triggered by a large-x proton would have 
correlated partners distinguishable from the background. In that respect the hadronization 
problem is similar to that at intermediate px in heavy-ion collision at LHC, where so many 
semi-hard jets are produced that shower partons are dense and can recombine with large 
p/tt ratio [H]. For the same reason as at large x studied here, it was also predicted that for 
triggers in the 10 < < 20 GeV/c range no correlation structure of associated particles 
would be found. Thus to a certain extent what we can learn about forward production at 
RHIC may reveal some aspects of the characteristics of what may be observed at intermediate 
Pt at midrapidity at LHC. 
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Figure 1: The x distributions of produced proton, pion and antiproton in Au-Au collisions 
at 6 = 3.3 fm and = 62.4 GeV for k = 0.7. 
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Figure 2: A plot of antiproton/proton ratio versus proton/pion ratio for various fixed values 
of X (solid lines) and k (dashed lines) for Au-Au collisions at 0-10% centrality. The theoretical 
curves are determined by calculating Hhr{x, k)/ Hh{x, k). The experimental point is from the 
BRAHMS data on Rpjp [S] and Rpj^ at = 62.4 GeV and ?7 = 3.2 ± 0.2. 
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Figure 3: The pt dependence of proton-to-pion ratio in Au-Au collisions at rj = 3.2. The 
dashed line is obtained from the ratio Hp[x{pT)]/ H.,^[x{pt)] with k, = 0.67. The solid line 
includes the factor V^(pt)/K-(pt)- The data (preliminary) are from [Sj. 



15 



I 



> 
o 



"D 

h- 
Q. 
"D 



"D 



10^ 



10^ 



-1 



-10 

Q. 

^ 10 



10 



-3 





r|=3.2 




0-10% 


■ ■ BRAHMS 




(h++h")/2 




■ — K=0.67 





0.5 



1 1.5 

(GeV/c) 



2.5 



Figure 4: The distribution of average charged hadron in Au-Au colhsions at = 3.2. The 
data are from [7j. The sohd hne is obtained by use of Eq. for /i^ = [p+p+1.2(7r++7r^)]/2 
and K = 0.67. 
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Figure 5: The pt dependence of antiproton-to-proton ratio for 0-10% centrality in Au-Au 
collisions at = 62.4 GeV, 77 = 3.2 and k = 0.67. 
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